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O Abstract 

^ jWe comment on the recent results of Refs. 01 01 on the bipolaron problem derived using an approximation of Gross - Tulub. It is 
Qh proved that, contrary to the claim made in Refs. [00, the bipolaron ground state energy calculated there in the strong- coupling 
<^ approximation has not been shown to constitute a variational upper bound. 
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1. Introduction 

Frohlich-bipolarons play a role in the study, e.g., of elec- 
tronic and magnetic properties of polar solids [0 §] |(1 
ana (nun 11311 . Bipolarons have, e. g ., been invoked in stud- 
ies of high- T c superconductivity JHO. It was shown in 00] 
(see also Ref. yj]) that bipolarons can exist only in a highly 
restricted stability -region of the (U,a) plane (where a is the 
electron-phonon coupling constant and U is the strength of the 
, Coulomb repulsion) and that some of the high-!T c oxides belong 
to this stability region, whereas conventional polaron-materials 
(alkali-halides, Silver-halides etc. . . ) have {U, a) parameters far 
away from the bipolaron stability region. Bipolarons are only 
stable at sufficiently large a, so that intermediate and strong 
coupling are relevant for the present discussion. Importantly 
for the study of high-r c superconductivity, it was found (a |7[] 
that a bipolaron binds more easily in 2D than in 3D and that 
the average pair-radius is a few Angstroms. The groundstate 
properties of bipolarons were studied further, e. g., in Refs. 

Recently, the large-bipolaron problem was approached in 
Refs. |0, 0] using the Gross - Tulub (GT) approximation 
iffl [I5I1 in the strong-coupling limit. Surprisingly low upper 
bounds to the bipolaron groundstate energy were arrived at in 
|0 0. In the limiting case for the ratio of the dielectric con- 
stants 77 = fioo/eo = and the Frohlich coupling constant a » 1, 
the bipolaron groundstate energy proposed as an upper bound in 
|0] is E bip (a » 1,77 = 0) ~ -0.414125a 2 , significantly lower 
than any other result in the literature. 

In the present communication we analyse the method and the 
results of Refs. (00 and we demonstrate that the approxima- 
tion for the bipolaron strong coupling groundstate energy ar- 
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rived at in (0 has not been shown - contrary to the claim in (0 
- to constitute an upper bound. 

2. General treatment in Refs. 100] 

Consider the two-polaron (bipolaron) system with the Hamil- 
tonian 
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with the coupling parameters 
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Here, a is the Frohlich polaron coupling constant, is the 
high-frequency dielectric constant. The system of units is cho- 
sen with 2m = 1 , h = 1 , and the LO-phonon frequency oiq = 1 . 
After the transformation to the center-of-mass and relative co- 
ordinates 

R= — 2 — > r = r i- r 2. (3) 
a first Lee-Low-Pines (LLP)-type canonical transformation 
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and the averaging of the transformed Hamiltonian with a trial 
wave function ip(r) for the relative motion (in the oscillatory 
form as in Refs. |00), 
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we arrive at the reduced Hamiltonian 



H 
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The second canonical transformation 



(6) 



(7) 



with the trial phonon shifts f k chosen as in Ref. J2l with the 
variational parameter p, 
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and the generalized Bogoliubov transformation used in Ref. 
| I5J result in the following variational functional for the bipo- 
laron groundstate energy 
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with the parameters 



a — 4 V2a, a — 
and the recoil energy 
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E R = —a 2 [l + Q(a,a)]. 
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which results from the recoil term |16| ^(^k^a^k) i n me 
Hamiltonian ©. The function 2 (a, a) is given by the integral 
expression IU5I1 



Q (a, a) 
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with the functions 
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and the parameters 
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3. Using the approximation of Ref. OH] in Ref. |0] 

In Ref. JH], the function g (a, a) [denoted in 1 15| as q(l/A)] 
has been replaced by the value Q m as 5.75. If the recoil energy 
<TTT~b is calculated using Q (a, a) = Qoo, one arrives at the result 
of Ref. H for the bipolaron energy (Eq. (15) of Ref. |H). 
However, this approximation does not guarantee a variational 
upper bound for the polaron and bipolaron groundstate energies 
for the following reason. 




Figure 1 : Q (a, a) given by Eq. 1121 



For finite parameters a and a, Q (a, a) is an increasing func- 
tion of both a and a, as shown in Fig. 1. Also in the strong- 
coupling limit, when both a and a tend to infinity, the asymp- 
totic expression for this function, 
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increases monotonically. The value ~ 5.75 has been ob- 
tained in Ref. 11511 assuming a finite cutoff for the phonon wave 
vectors. Indeed, the phonon wave vectors are restricted by a 
short-wavelength cutoff at the boundary of the Brillouin zone. 
Consequently, in the strong-coupling limit a — » oo, the posi- 
tion of the steep maximum mentioned in Ref. (15] Q can lie 
beyond the integration range. However, as far as the calcula- 
tion of the (bi)polaron energy is performed within the contin- 
uum approach, it is only consistent either to avoid a cutoff all 
together (in the continuum formalism), or to introduce a cut- 
off from the very beginning of the calculation. A consistent 
treatment of the polaron problem in the Gross - Tulub approxi- 
mation using a cutoff was performed in Ref. lfl6ll . It was found 
that a cutoff leads to the appearance of additional positive terms 
in the polaron recoil energy. These terms were not found in 
Ref. 1 1511 . Being positive, they definitely lead to an increase of 
the groundstate energy. Therefore missing these terms can lead 



1 As written in Ref. 11511 on p. 4, "It is of interest to note that as A — » oo the 
integrand in (2.12) has a steep maximum at y 4 = 3/1/4; however, if we take into 
account that the domain of integration over y is in fact limited and if we use the 
values g 2 a 10 considered in the following, this singularity does not arise." 
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to the violation of the variational principle. Hence the bipo- 
laron groundstate energy arrived at in Refs. [0 0] is incorrectly 
claimed to constitute an upper bound for the groundstate energy 
of the bipolaron. 



4. Bipolaron variational energy in the continuum approach 

Let us consistently consider the polaron and bipolaron 
groundstate energies within the continuum approximation us- 
ing the complete recoil energy (i. e, without a cutoff). In the 
strong-coupling limit, using $9^ with the asymptotic expression 
( TTol ). the polaron and bipolaron groundstate energies are then 
found to be 
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Therefore, when the (bi)polaron groundstate energy is consis- 
tently calculated within the continuum GT approach, an incor- 
rect dependence E oc or' 3 results in the strong-coupling limit. 
This problem, however, was not realised by the author of Refs. 

Hi. 
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Figure 2: Bipolaron variational groundstate energy calculated using the com- 
plete recoil energy within the Gross - Tulub scheme in the continuum approach. 
For comparison, the bipolaron groundstate energy calculated in Ref. (g| and 
twice the groundstate energy of a single polaron from Ref. fnll are plotted in 
the same graph. 



In order to check whether the continuum variational GT 
approach is adequate for intermediate a, we have calculated 
the bipolaron groundstate energy using the expression (|9) with 
Q (a, a)given by <TT~2T > instead of the value Q M « 5.75 used in 
Refs. [1LI2U- In Fig. 2, this bipolaron groundstate energy is 
compared with the variational result of Ref. H and with twice 
the energy of a single polaron calculated using the path-integral 
variational method of Feynman Jl7ll . It is seen that for interme- 
diate a, the bipolaron groundstate energy calculated using the 
variational GT scheme in the continuum approach lies above 
twice the energy of a single polaron. 



5. Conclusion 

In the present communication we have proved that the strong- 
coupling expression for the bipolaron groundstate energy c alcu- 
lated in Refs. contrary to what is claimed in those works 

is not justified as upper bound for the bipolaron groundstate en- 
ergy. A contribution to the recoil energy due to a momentum 
cutoff Jl6Tl within the GT scheme is missed in Q],Q]. Thus the 
results of Refs. 01 0] have been obtained using the incomplete 
recoil energy. 

We have also worked out the GT scheme using the complete 
recoil energy within the continuum approach as is necessary 
for a consistent theory. This leads to an incorrect dependence 
Ebip c 4 ^ 3 (instead of Ebip oc a 2 ) for the polaron groundstate 
energy in the strong-coupling limit. A consistent calculation 
within the continuum approach, however, was not performed in 
Refs. Hi]. 

Note also that for intermediate a, the continuum GT method 
leads to a bipolaron groundstate energy higher than twice the 
groundstate energy for a single polaron. Therefore this method 
fails to describe a bound bipolaron state for any a. 
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